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Abstract 

Hamiltonians of a wide-spread class of strongly coupled quantum system models are 
r*"""- . expressed as nonlinear functions of sl(2) generators. It enables us to use the si (2) for- 

malism, in particular, sl(2) generalized coherent states (GCS) for solving both spectral 
and evolution tasks. In such a manner, using standard variational schemes with sl(2) 



GCS as trial functions we find new analytical expressions for energy spectra and non- 
linear evolution equations for cluster dynamics variables in mean-field approximations 
which are beyond quasi-harmonic ones obtained earlier. General results are illustrated 
on certain concrete models of quantum optics and laser physics. 
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2 : 1 Introduction 



For last decades a great attention has been paid to solve and to examine different dynami- 
cal problems for quantum strongly coupled systems whose interaction Hamiltonians are ex- 
pressed by nonlinear functions of operators describing subsystems (see, e.g., [1-9] and refer- 
Ch ences therein). However, as a rule, for these purposes numerical calculations are mainly used 
^ ■ while analytical techniques available either deal with special forms of model Hamiltonians 
O" 1 ' (including their different semiclassical versions) and initial quantum states [1-5,7-9] or require 
>• . lengthy and tedious calculations (as it is the case, e.g., for the algebraic Bethe ansatz [6]). 

Recently, a new universal Lie-algebraic approach has been developed [10-12] to get solu- 
^ [ tions of both spectral and evolution problems for some nonlinear quantum models of strongly 
coupled subsystems having symmetry groups Gi nv . It was based on exploiting a formalism of 
polynomial Lie algebras g p d as dynamic symmetry algebras g DS of models under study with 
generators of these algebras g p d being Gj n „-invariant collective (cluster) dynamic variables in 
whose terms model dynamics are described completely. (In fact, such a reformulation of orig- 
inal problems in terms of ^-variables is similiar to the well-known procedure of exclusion of 
cyclic variables in classical mechanics [13].) Specifically, this approach enabled us to develop 
some efficient techniques for solving physical tasks in the case of g DS = sl p d(2), when model 
Hamiltonians H are expressed as follows 

H = aV + gV + + g*V- + C, [V a ,C] = 0, V^ = (V + ) + , (1.1) 

where C = C([Ri\) is a function of a set of commuting operators (model integrals of motion) 
Ri, i — 1, 2, ... and Vo, V± are the sl p d(2) generators satisfying the commutation relations 

[Vo, V ± ] = ±V ± , [VI, V+] = (V + 1) - *(V ), *(Vo) = A J] % + H{Rj})), (L2a) 

i=l 
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[*(Ro),V a ]=0 Va = 0,±, *(Ro) = *(V )-V + V-, (1.2b) 

where n(\l/) is the polynomial ^ degree in the variable Vo, ^{Rq) is the sl p d{2) Casimir 
operator (with R being the "lowest weight operator") and hereafter the identity operator 
symbol I is omitted in expressions like ^{V + al). The structure polynomials ^(V ) depend 
additionally on {Ri,i — 1, . . .}, and their exact expressions for some wide-spread classes of 
concrete models were given in [10-12] (see also Section 5). 

All techniques [10-12] essentially use expansions of evolution operators Un{t), generalized 
coherent states (GCS), energy eigenfunctions \Ef > and other important physical quantities 
by power series in the sl p d{2) shift generators V± as well as commutation relations (1.2) and 
the characteristic equation 

(V+V- - *(V ) = -*(Ro))\l(h) = (1.3) 

fulfilled on Hilbert spaces L(H) of quantum model states due to the complementarity of 
groups G inv and algebras sl pd {2) [10]. Specifically, Eq. (1.3) implies a spectral decomposition 

L{H) = £ L([k}), L([l l ])^Span{\[l l ];v)=Af([k};v)V;\[l t ]}}, 

Pi] 

V \U,v) = (l + v)\[h};v), R i \[l i ];v) = l i \[l i ];v),i = 0,l,..., V_\[h}) =0 (1.4) 

of spaces L(H) in direct sums of the subspaces L ([/j]) = L(l ,li,...) which are irreducible 
with respect to joint actions of algebras sl p d{2) and symmetry groups G ills ; lowest weights Iq 
depend on other quantum numbers — 1, . . . due to the relation \l/(/o) = implied by Eq. 
(1.3). From the physical point of view, the decomposition (1.4) means that the model Hamil- 
tonian matrices in the symmetry adapted orthonormalized bases { | [/j] ; v) } have block-diagonal 
forms and subspaces £([/«]) describe specific "s/ p d(2)-domains" evolving independently in time 
under action of Hamiltonians (1.1). We also distinguish compact (sw p d(2)) and non-compact 
(su p d(l, 1)) versions of sl p d(2) algebras depending on whether dimensions d([li\) of the spaces 
L([/j]) are finite or infinite. 

Then, using restrictions H^.j of Eq. (1.1) on £([/«]), simple algebraic calculation schemes 
were developed for finding evolution operators U H [t) = S/L- C o^ / + u f(V ;t), amplitudes 
Q v (Ef) = ([li];v\Ef) of energy eigenstates \Ef) expansions in orthonormalized bases {\[h}; v)}) 
and appropriate energy spectra {Ef} of bound states [10]. In the paper [11] some explicit 
integral representations were found for amplitudes Q V (E), eigenenergies {E a } and "evolution 
coefficients" Uf(V ;t) with the help of a specific "dressing" (mapping) of solutions of some 
auxiliary exactly solvable tasks with the dynamic algebra s/(2). 

However, all exact results obtained do not yield simple working formulas for analysing 
models (1.1) and revealing different physical effects (e.g., a structure of collapses and revivals 
of the Rabi oscillations [2,8], bifurcations and singularities of quasiclassical solutions [5] etc.) 
at arbitrary initial quantum states of models. Therefore, it is necessary to develop some 
simple techniques, in particular, to get some closed, perhaps, approximate expressions for 
evolution operators, energy eigenvalues and wave eigenfunctions, which would describe main 
physical peculiarities of model dynamics with a good accuracy (cf. [5,8,9]). Below we examine 
some possibilities along these lines for models (1.1)-(1.4) by means of reformulating them in 
terms of the usual sl(2) algebra formalism and developing variational schemes corresponding 
to quasiclassical approximations (QAs) for these models by analogy with developments [5,14- 
16]. 
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The work is organized as follows. In Section 2 we first reformulate models (1.1)-(1.4) in 
terms of the usual Lie algebra si (2) formalism, and then discuss possibilities of extending 
the standard s/(2)-techniques to analyse such reformulated models. In Section 3 a scheme 
is given for obtaining QAs of these models by using variational principles [5,17] and energy 
functionals constructed with the help of the SL{2) group GCS [16]; these QAs are new for 
original models because they take into account a strong coupling of interacting subsystems 
in contrast with standard QAs. Specifically, in such a manner new analytical expressions 
are obtained for energy spectra which are essentially non-equidistant on each subspace £([/;]) 
with its dimension d([li\) > 4. In Section 4 we discuss such approximations for a quasiclassical 
description of dynamics of s/(2)-clusters (characteristic model exicitations) and time evolution 
of uncoupled dynamical variables; specifically, nonlinear evolution equations of the Bloch type 
are obtained for s/(2)-cluster variables. In Section 5 a specification of general results is given 
for a class of models widely used in quantum optics and laser physics. In conclusion some of 
prospects of developing this approach are discussed. 

2 A nonlinear sl(2) formulation and a general opera- 
tor analysis of quantum models with linear sl p d(2) dy- 
namic algebras 

We can reformulate models (1.1), (1.4) in terms of si (2) generators using a realization of 
the sl p d{2) algebras in terms of special elements of extended enveloping algebras Uy(sl(2)) of 
the familiar algebra si (2) [12]. This realization is established via the generalized Holstein- 
Primakoff mapping [10] 

Y = V -B«Tj, Y + = V + IHY )]-^, $(Y ) = J^ff^tvV Y ~ = ^ 

(JtY )(±J+ 1 + Y ) 

[Y , Y ± ] = ±Y ± , [y_, Y + ] = T 2Y (2.16) 
where Y a are the si (2) generators, R , q=J are lowest weight operators of the sl pd {2) and 
si (2) algebras respectively: ^(Rq))\l(h) = 0, [J(J ± 1) =F Y+Y- ~ Yq^]\l(h) = 0. Appropriate 
specifications of Eqs. (2.1a) on subspaces £([/»]) are obtained by the substitution R Q — > l Q , J — > 
J and hereafter upper/lower signs in (2.1) corresponding to the su(2)/su(l, 1) algebras are 
chosen for finite/infinite dimensions <i([/i]) of the spaces L([/j]). Note that, by definition (2.1), 
functions $(lo) on L([li\) are polynomials of the (n(^) — 2)-th degree in the variable Y at 
relevant values of J [10,12]. 

Then, using Eqs. (2.1) one may re-write Hamiltonians (1.1) in terms of Y a as follows, 

H = aY + Y + g(Y ) + g + (Y )Y_ + C , 

g(Y ) = gjMYo), C = C'{[Rl J) = C([R t ]) + a(R ± J) (2.2) 

Restrictions Hua = PuaH of Hamiltonians (1.1) on spaces L ([/j]) (with Pry = J2 V \[h]] v )( v 'i [h]\ 
being appropriate central projectors) are obtained by the substitution Rq — > Iq, J — > J in 
Eq.(2.2). Respectively, basis vectors |[/j];f) of spaces £([/»]) are given in terms of Y a as 
follows, 

\[l i ];v)=tf(J,v)(Y + y\[l i ]) (2.3) 
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where M~ 2 (J,v) = v\(2J)\/(2J-v)\ for su(2) and Af~ 2 (J,v) = v\T(2J +v) /Y(2J) for su(l, 1). 
Evidently, Eq. (2.2) resembles Hamiltonians of semi-classical si (2) "linearized" versions of 
matter-radiation interaction models [4,8,9,12] but with operator (intensity-dependent) cou- 
pling coefficients g(Y ) (cf. [3,4,7]). Emphasize, however, a collective (not associated with a 
single subsystem) nature of operators Y a that leads, when substituting g(Yo) in Eq. (2.2) by 
an "effective coupling constant", to a non-standard ("cluster") QA of original models [10] dis- 
tinguished from standard semi-classical limits [8,18,19] where a part of interacting subsystems 
is described classically. 

If n(^) = 2, then $(Y ) = l,sl p d(2) = sl(2),R = =fJ, and we have a powerful tool for 
solving both spectral and evolution tasks yielded by the GCS formalism [16] related to the 
SL(2) group displacement operators 

S Y (£ = re w ) = exp(£F + - f*y_) = exp[t(r)e id Y + ] exp[-2 In c(r)Y ] exp[-t(r)e _ifl y_] = 

£ Y{SJ(Y ;0, Y- k = Y*([* 2 (Yo)] W ) for k > (2.4a) 

f=-oo 

where t(r) = tanr/ tanhr, c(r) = cosr/ coshr, s(r) = sinr/ sinhr for su(2)/su(l, 1) and 
[^ 2 (y )] (fc) = (±l) fe (±J + F ) (fc) (±^ ~Y + k) {k \ = A(A - 1)...(A -x + 1) (2.46) 

Sj(Y 0] = 1 )] 2F 1 {TJ ~ Y , -Y ± J + 1; / + 1; ±[s(r)] 2 ) exp[-2 In c(r)Y ] (2.4c) 

with 2 -P\(---) being the Gauss hypergeometric function [20]. 

Specifically, in this case, using the well-known sl(2) transformation properties of operators 
Y a under the action of [16,12], 

S Y (0Y + S Y (^ = F + (0 = [c(r)] 2 Y + ± e^[s(2r)Y - e -*[ a (r)] 2 y_], 

H)WO f = ^o(0 = c(2r)Y - 6 -^[e* e Y + + e~ ie Y_, Y.® = (Y + (0)\ (2-5) 
Hamiltonians H can be transformed into the form 

H(0 = S Y (0HS Y (^ = C' + Y A (a, g; £) + Y + A + (a, g; £) + Y_A* + (a, g; £) (2.6a) 

At the values £ = of the parameter £ with tan2r/ tanh2r = ^ for su(2)/su(l, 1) one 
gets A + (a,g;£) = 0, and the Hamiltonian #[y(£) takes the form 

#(6) = C" + r 0v / « 2 ±4|^| 2 (2.66) 

which is diagonal on eigenf unctions |[Zj];i>). Therefore, original Hamiltonians H have within 
each £([/i]) equidistant spectra with eigenenergies 

E([l i ];v) = C + (TJ + v) y Ja>±4\g\\ C = C ([% J) (2.7a) 

and eigenfunctions 

\[li];v;Zo) = S Y (^\[k];v) =eM-ZoY + + e o Y-)\[k];v;) = £S / ,(J;ff ) r)|[I i ];/) ) 
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S fv (J; 9, r) = (f-J)W( JJ) 2Fli ~ V > ± 2 J + 1; / - t; + 1; ±[s(r)}% 

(2-76) 

where t(r) = ± (—a + \Ja 2 ± 4|g| 2 ) /2|c/| and A/"(J, ...) are normalization constants from Eq. 
(2-3). 

Similarly, if sl p d(2) = sl(2), operators Sy (£(£)) are "principal" parts in the evolution 
operators Un(t) = exp(ia(t)Yo)SY{£{t)) with a (£),£(£) being c-number functions in t which 
are determined from disentangling the exponent exp(^iif) (or, when g are time-dependent 
functions, from a set of non-linear differential equations corresponding to classical motions) 
[16,19]. 

However, for arbitrary degrees n of polynomials \I/(Vb) Hamiltonians (2.2) are essentially 
nonlinear in si (2) generators Y a , and, therefore, the situation is very changed. Specifically, in 
general cases it is unlikely to diagonalize H with the help of operators <Sy(0 since analogs of 
Eq. (2.6a), 

H(0 = S Y (OHS Y (^ = aY (O + Y + (t)g(y (t)) + g + (Y (£))Y_(£) + C,' (2.8) 

and even their restrictions H[ii\(£) = S , y(^)f/"[; i ]S'y(^) 1 ' on multi-dimensional spaces L([li\) 
contain (after expanding g(Yo(£)) in power series) many terms with higher powers of Y±. The 
task is also not simplified when using in Eq. (2.8) operators <Sy(£) = exp(£V+ — instead 
of <Sy(£) because we have not suitable analogs of the "disentangling theorem" (2.4) and finite- 
dimensional transformations (2.5) for operators 5V(0 [12]. Therefore it is necessary to use in 
Eq. (2.8) more general (perhaps, non- or multi-parametric) forms of diagonalizing operators 
S given, e.g., by power series 

S= E Y +Sf(Y ), Y- k = Y k ([^ 2 (Y )]^Y VA;>0 (2.9) 

f=-oo 

with undetermined (unlike Eq. (2.4)) coefficients Sf(Y ) and satisfying the unitarity condi- 
tions S& = S^S = I. 

Substituting Eq. (2.9) in the scheme (2.8) one gets after some algebra nonlinear analogs 
ofEqs. (2.6a) 

OO 

H = SH& = C'+ E y+hf(Y ) (2.10a) 

f=—oo 

where 

OO 

hf(Y )= E ^ 2 (Y )]^S* k (Y -k)[a(Y -k)S k+f (Y -k) + 

k=—oo 

gy/$(Y - k)S k _ 1+f (Y -k + l) + g*y/$(Y -k-l)S k+1+f (Y -k- l)* 2 (y - k)} (2.106) 

and h_ f {Y Q ) = h}(Yo - /)[* 2 (y )] (/) , [^(^o)]^ = ([* 2 (^o + /)] {/) )^ for all / > 0. 
The conditions hf(Y ) = for all / ^ yield nonlinear analogs of Eqs. (2.6b), (2.7b), 

OO 

H° = S°HS°* = C' + h° (Y ) = J2 [MYo)] {k) S° k *(Y - k)[a(Y - k)S° k (Y - k) + 

fc = — OO 

g^iYo-kjS^Yo -k + l)+ g*y/$(Y -k-l)S k+1 (Y - k - 1)* 2 (Y - k)], (2.11a) 
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S([Z i ]; U )=C+([Z i ]; U |/ l 8(y )|[ii];u) = C + / l g(iFJ + u) (2.115) 

expressed in terms of the coefficients Sf(Y ) which satisfy Eqs. (2.10b) with hf(Y ) = for 
/ 7^ and simultaneuosly are solutions of the set of algebraic operator equations, 

l-aY + h (Y + f)]S° f (Y ) = 

gJ^Yojs^Yo + 1) + g*y/$(Y - l)S° f+1 (Y - 1)* 2 (Y ) V|/|>1 (2.12) 

resulting from the condition S°H = [C + h%(Y )]S° [12] where S° is given by Eq. (2.9) with 
"coefficients" S° f (Y ). 

In the case of $(lo) — 1 Eqs. (2.12) are solved in terms of hypergeometric functions 
2-Fi(...) as it follows from Eq. (2.4c), but in general they, probably, determine certain g-special 
functions due to relations of sl p d(2) algebras with certain g-deformed algebras [21]. Without 
using any specifications of operators S, due to the relation Q v (Ef) = Sj_ v (±J + v) jffj^ , the 
task of solving these equations is equivalent to that for finding amplitudes Q v {Ef) related 
to new classes of orthogonal functioms [10]. Note that this task is simplified in the compact 
(su(2)) case, when all subspaces £(&]) have finite dimensions cZ([Zj]) = 2J + 1, and all series 
in Eqs. (2.9)-(2.11) are terminating due to Eq. (1.3) and the relation (y±) 2J+1 |L([y) — 0. 
Therefore, eigenf unctions \E V ) = S^\[li};v) may be represented by polynomials 

2J 

\K) = EQJXflfc]) = A V H(Y + - <)[[/,]), (2.13) 

f=0 r 

where amplitudes Q V f are expressed as symmetric functions in variables K v r : 

2J 

Q2J = A v , Q2J-1 = ~A V ^ K v r , 

r=0 

Q 2 J-f — ( — ^- Y A v ^ K r 1 K r 2 --- K r f ^ —■> Ql = ( — 1 ) 2 J ^ K ri K r 2 • • - K r 2 j ' (2-14) 
l<r 1 <r 2 <...<rf<2J 

At the same time eigenenergies E([li];v) determined by the boundary condition [10] 

[(Z + 2J)a - E([k];v) + C([h})]Q v 2J + gQ u 2J _ x = 0, (2.15) 
can be written down in the form 

2J 

E([h];v) =C+Ja-gJ2<, C = C([Z i ]) + (Z + J)a (2.16) 

r=0 

of a sum of 2 J + 1 spectral functions as it is prescribed by the algebraic Bethe ansatz [6]. 
In fact, Eqs. (2.13)-(2.16) give for models given by Eqs. (1.1)-(1.4) a new, G^-invariant 
formulation of this ansatz in terms of the su(2) algebra which is simpler and more efficient 
in comparison with its initial non-invariant version [6] because the algorithm [10] for finding 
amplitudes Q V f and eigenenergies E([li];v) does not require a preliminary determination of 
parameters k v t . We also note that in the resonance case (when a = in (2.2)), using Eqs. 
(2.12), one can get analytical solutions for amplitudes Q v j- in the form of multiple sums which, 
however, are not suitable for practical purposes. 
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So, direct generalizations of "linear" schemes (2.6) to the case of non-linear Hamiltonians 
(2.2) do not yield simple analytical formulas for exact solutions of spectral tasks; a similar 
situation is also with respect to evolution problems. Nevertheless, the formalism of the SL(2) 
GCS |[Jj];i>;£) = «Sy(£)t| [k];v) can be an efficient tool for analysing such models [5,10,14-16] 
and for getting approximate analytical solutions of both spectral and evolution problems. 
Specifically, a simplest example of such approximations was given in [10] by mapping (with 
the help of the change V a — > Y a ) Hamiltonians (1.1) into Hamiltonians H s ^ 2 ) which are linear 
in sl(2) generators Y a (but with modified constants a,g) and have on each fixed subspace 
L([/j]) equidistant energy spectra given by Eq. (2.7a). However, this (quasi)equidistant ap- 
proximation, in fact, corresponding to a substitution of certain effective coupling constants g 
instead of true operator entities g(Yo) in Eq. (2.2), does not enable to display many peculiar- 
ities of models (1.1) related to essentially non-equidistant parts of their spectra. Therefore, 
it is needed in corrections, e.g., with the help of iterative schemes [8,14,15]; specifically, one 
may develop perturbative schemes by using expansions of operator entities g(Y ) in Taylor 
series in Yq as it was made implicitly for the Dicke model in [8,9]. But there exists a more ef- 
fective, incorporating many peculiarities of models (1.1), way to amend the quasi-equidistant 
approximation. 



3 SL(2) quasiclassical approximations: energy function- 
al and variational energy spectra 

This way is in applying SL{2) GCS [[/«]; t>;£) from Eq. (2.7b) as trial functions in the varia- 
tional schemes [17] of determining energy spectra and quasiclassical dynamics [5,15]. Indeed, 
because of the isomorhism of quantum and quasiclassical dynamics for sl(2) linear Hamil- 
tonians [15,16], the results (2.7) can be obtained with the help of the variational scheme 
determined by the stationarity conditions 

for the energy functional H([li]',v,g) = ([li]',v;Z\H\[li];v',Z) = ([h]; v\C + Y A (a, g; f) + 
Y + A + (a, g;£) + Y_A* + (a, g; £)\[k]; v ) (cf. (2.6a)). Similarly, following the standard varia- 
tional approach [17,5], the calculation schemes (3.1) may be extended to the case of nonlinear 
Hamiltonians (2.2) by using the energy functional 

n c \U;v-i) = {[li];v;t\H\[li];v;t) = (v; [l t }\H(0\[k];v) (3.2) 

where superscript cq denotes "cluster" (strongly correlated) QAs (as contrasted with standard 
QAs dealing with weakly or non-correlated subsystems) and H(£) are given by Eq. (2.8) or 
Eqs. (2.10) with S f (Y ) = Sj{Y ;£) from Eqs. (2.4c). 

Note that, due to the Hermitian conjugacy relations Y + g(Yo) = (g + (Yq)Y-)\ the condition 
(3.1a) gives e ld = g/\g\ (as in the linear case); in fact, this condition (3.1a) can be eliminated 
at once by the simple gauge transformation Y a — > exp(— ieey)Y a , g = \g\ exp(ry), in Eqs. (2.2) 
which preserves commutation relations (2.1b). Furthermore, due to the form of trial functions 
and the unitarity of operators it is sufficiently to solve Eq. (3.1b) only for finding 

states \ [k};v = 0; £) and to use the only real root of one of Eqs. (3.1b) for all v that ensures 
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automatically the orthogonality of eigenfunctions. Naturally, results thus obtained are not 
expected to coincide with exact solutions on all subspaces due to an essential non- 

linearity of Hamiltonians (2.2b) and their non-equivalence (unlike Eq. (2.6b)) to diagonal 
parts of Eq. (2.10a); however, they yield "smooth" (analytical) solutions which are in a sense 
most close to exact ones (cf. [5,14]). Without discussing all aspects of such extensions we 
give below two approximations for energy spectra obtained by inserting in Eq. (3.2) H(£) 
given by Eqs. (2.10) and Eq. (2.8) respectively. 

In the first case, using Eq. (2.5) for Y (£ ), Eq. (2.7b) for and defining relations 

for the si (2) algebra, one gets the following "cluster" QAs 



E c «([k];v) = C + a(vT J)c(2r) + ^([Z i ]; U ;eomV $ ( y o)IM;u;6)} 



C + a(v =F J)c(2r) - 2\g\ ]T \S fv (J; g, r)S f+lv (J; g, r) \ yjv(l + 1 + /), (3.3a) 

/>o 

E^U 0) = C T a. M 2r) - 2\g\ [<r )]*» £ ^^J^) (3-36) 

for energy eigenvalues E c<1 ({k}; v) = H cq ({li];v;£ ) where C = C([k\) + a(l ± J), & = 
rg/\g\,2®{A} = A* + A, c(r) = cos r/ cosh r, ^(Z + 1 + /) = $( T J + f)(2 J T /)(/ + 1), 
Af(J,...) are normalization constants from Eq. (2.3) and functions Sf v (J;g,r) are given by 
Eq. (2.7b) but with values of the parameter r determined by real solutions of the algebraic 
equation 



a(1 ± " } = ^ mrwwTT) [TAa * + ± ° )< 2 ' + ^ 

a = t(r) = tanr/ tanhr (3.4) 

which follows from Eqs. (3.1b) and (3.3b). Obviously, unlike the linear case, diagonaliz- 
ing values of r depend on both constants g, a and quantum numbers U labeling G-invariant 
subspaces L ([/«]). 

In the second case it is difficult to obtain exact analytical formulas like Eq. (3.3) due to 
presence of square roots in Eq. (2.8). However we can get another approximation for energy 
spectra if replacing the energy functionals (3.2) by their (corresponding to the Ehrenfest 
theorem with respect to cluster variables Yj) mean-field approximations 

n cmf ([k] ; v; = a < F (0 > +2K{< Y + (0 > g(< F (0 >)} + C, 

<Y a (0>=(v;[k]\Y a (0\[k];v), 2^{A} = A + A* (3.5) 

Then, inserting Eqs. (2.5) in Eq.(3.5), one finds the "cluster" mean-field approximations 
E cm * ( [I j] ; v) for eigenenergies, 

E cmf ([% v) = C + a(v =f J)c(2r) - 2\g\(J t v)s(2r) ^$((=F J + v)c(2r)) (3.6a) 

where r is determined from the equation 

o , n , , ^ I , T ~ ~ J[s(2r)] 2 $'(^Jc(2r)) T /o xx d$(ar) . 

— s 2r = ±2c(2r) v /$( T Jc(2r)) + [ 1 >\ y J $' ( T Jc 2r = -^L= T j c (2r) 

(3.66) 
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Let us make some remarks concerning results obtained. 

Remark 1. As is seen from Eq. (3.3), its general structure coincides with the energy for- 
mula given by Eq. (2.16), and spectral functions \Sf v (J; g, r)Sf + i v (J; g, r)\yJ^(lo + 1 + /) = 
Ef(r; J; v) are nonlinear in the discrete variable v labeling energy levels within that P ro " 

vides a non-equdistant character of energy spectra within fixed subspaces £([/;]) at d([li\) > 3. 
Besides, due to the presence of square roots in Eqs. (3.3), (3.6a) different eigenfrequencies 
u) v = E([li];v)/h are incommensurable: mu Vl ^ nu V2 that is an indicator of complex dy- 
namics manifesting in such phenomena as collapses-revivals of the Rabi oscillations [2,8] and 
singular and pre-chaotic dynamic regimes in phase spaces of models [13,22,23]. Evidently, 
it is hardly possible to obtain these features of models by using GCS related to uncoupled 
subsystems (cf. [4,8,18] and Section 5 of the present paper). 

Remark 2. In the compact (su(2)) case the r.h.s. of Eq. (3.4) is a polynomial of the degree 
2 J + 1, and, in general, Eq. (3.4) may have 2 J + 1 different roots corresponding to 2 J + 1 
different stationary values of the energy functional W([Zj]; v; £). Therefore, one may assume 
that it is possible to get more simple expressions for certain E([li];v) using £([/;]; 0) from 
Eq. (3.3b) with different real roots r« of Eq. (3.4); specifically, it is the case for dimensions 
rf([/j]) = 2 when Eqs. (3.3)-(3.4) give exact results. However, using the well-known expressions 
for overlap integrals of SU(2) GCS [16], one can show that in general only two £77(2) GCS 
with different real roots may be mutually orthogonal. 

Remark 3. Obviously, Eq. (3.3) generalizes Eq. (2.7a) for the (quasi)equidistant approx- 
imation abovementioned. Indeed, when replacing the functions ^(^fJ + /) by their certain 
(and the same for all labels v) "average" values, series in (3.3), (3.4) are summed up, and 
Eq. (3.3) is reduced to Eq. (2.7a); Taylor series expansions of functions \J $(^7 + /) provide 
perturbative corrections related to higher degrees of the an-harmonicity of Hamiltonians (2.2). 
At the same time Eqs. (3.5)-(3.6) yield an intermediate (related to a more fine "averaging" 
procedure (3.5)) approximation retaining the main characteristic feature of Eq. (3.3) (a non- 
equdistant character of energy spectra within fixed subspaces L ([/;])) but being simpler in its 
form that is important from the practical point of view. Besides, Eqs. (3.4) and (3.6b) are 
simplified in the resonance case when a = 0. 

Remark 4- In fact, solving Eqs. (3.1b) one can get a whole series of competitive potential 
solutions (corresponding to different roots r\ and v ) which may approximate exact ones with 
a good accuracy in particular parts of energy spectra. (This situation resembles that occuring 
in the stationary phase calculations of the path integral approach when one needs to take into 
account contributions of several classical trajectories [24,25].) A final selection of the most 
adequate value r may be made with the help of a "quality criterion" of QCAs on subspaces 
L([/j]). For example, one can estimate an accuracy of QCAs obtained by means of the "energy 
error" functionals [10] 

§ P i]{HjH c q /cm f) = |Tr W ( H - H <*/°»f mTr [h]{Hn p= 1,2 (3.7) 
giving "energy-trace" proximity meausures of the exact Hamiltonians (2.4) and their QCAs 
H^iiYt}-^) = ]T E*([U];v)\[li];v;t )(fav- t [h}\ = C + S Y (^ ) + h f (Y ; £ )Sy(£o), (3-8a) 

[h],v 

# cm/ (TO;6) = E ^ cm/ (M;^)IW;^;6>(6;^ [h}\ = 

[h],v 
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C + S Y (i«) + [a(Y Q )c(2r) ± 2|^|F oS (2r) [4> n _ 2 {{Y )c{2r))Y' 2 \S Y {io) (3.86) 

on subspaces L([Zj]); Tr uaA = Y^v( v i[k]\A\[li];v). Furthermore, functionals (3.7) may be 
used in alternative "minimization schemes" of determining the paprameter r$. 

4 Variational quasiclassical dynamics of 5L(2)-clusters 
and time evolution of uncoupled variables 

The energy functionals (3.2) and their mean-field approximations (3.5) may be also used for 
a quasiclassical analysis of time evolution of cluster dynamical variables related to the si (2) 
generators Y { (cf. [5]). As is known, when Hamiltonians (2.2) are linear in si (2) generators, 
quasiclassical dynamics is isomorphic to the exact quantum one [14-16] and is described by 
the classical Hamiltonian equations [5,14,16] 

q=^, P=-fp H = {z{t);[l i ]\H\[%z(t)) (4.1a) 

for "motion" of the canonical parameters p,q of the SL{2) GCS -z(i)) — exp(— z(t)Y+ + 
z(t)*Y-)\[li\) as trial functions in the time-dependent Hartree-Fock variational scheme [17] 
with the Lagrangian C = (z(t); [k]\(id/dt - H)\[k]; z(t)); q = 9, p = (z(t); [k] \Y \[k\- z(t)) = 
=FJc(2r), z = rexp(— i#). An equivalent formulation can be given in Y = (Yi,Y 2 ,Y ) space 
using sl{2) vector Euler-Lagrange equations [5], 

y = \VH x VC, y = {y 1 ,y 2 ,yo), Vi = (z(t); [k] |Y^]; z(t)), V = {d/d yi ,d/dy 2 ,d/dy ), 

C = ±yl + yl + y 2 2 = ±J 2 , y± = y 1 ±iy 2 , A x B = (A 2 B -A o B 2 ,A o B 1 -A 1 B ,A 1 B 2 -A 2 B 1 ) 

(4.16) 

which yield linear quasiclassical Bloch-type equations for sl(2) linear Hamiltonians [19]. 

In the general case of nonlinear Hamiltonians (2.2) Eqs. (4.1) with H given by Eqs. 
(3.2) and (3.5) at £ = z*,v = also describe a quasiclassical SL(2) "cluster" dynamics 
of models under study which, however, is not isomorphic to the exact quantum one [5]. 
Besides, Eqs. (4.1) obtained with the help of GCS |[/j];z(t)) describe dynamics of SL{2) 
clusters within each subspace L([U]) separately. Specifically, Eqs. (4.1a) determine a "sl(2) 
linearized" QA U l H .^(t) oc P[y exp(— z(t)Y + + z(t)*Y_) of L([/j])-restricted evolution operators 
UH;\k](t) = P\ii]UH(i) when initial wave functions \ipo) are equal to |[/j]) (cf. [19]); in a 
sense, this approximation is equivalent to that obtained by substitutions in Eq. (2.2b) time- 
dependent coupling functions g{t) (compatible with solutions of Eqs. (4.1a)) instead of g(Y ) 
(cf. [4,23]). However, for general initial wave functions |^ ) £ L(H) it is necessary to 
generalize these equations, e.g., by using GCS exp(— z(t)Y + + z(t)*Y^)\ip ). Without dwelling 
on a detailed analysis of this topic we write down examples of Eqs. (4.1) when appropriate 
explicit expressions for TC are obtained from Eqs. (3.3b) and (3.6a) at v = by means of the 
substitutions 



2\ g \ - ge -^ + g*e^ = 2[^{g} cosg + <3{g} sing], 2iZ{g} = g - g\ 2^{g} = g + g*, 

-Vi • V2 

i =, sm q = — = 



=pJc(2r) = y = p, cosg = ^ , sing = ^ 2 (4.2) 
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where the first line is taken from the substitution £ — r 9/\g\ ~^ z * — rexp(ig) in Eqs.(3.3) 
and the second one is a direct consequence of Eqs. (2.5). 

Then, from Eqs. (3.3b), (4.1b) and (4.2) one gets essentially nonlinear quasiclassical 
Bloch-type equations 

y = \vn cq x V C, SjH^ = (2^{g}e(y ), -2S{g}e(y ),a + 2[to{g}y 1 -S{g}y 2 ]de(y )/dy ), 

Qiyo) =[—) g mW7TW?g - VC = 2 (yi , y2 , ±yo ) (4.3) 

At the same time, using substitutions (4.2) in Eqs. (3.6a) and inserting them in Eqs. (4.1) one 
finds in the mean-field approximation (3.5), respectively, canonical Hamiltonians equations 

q = a T cosg + %{g} sin q] [±( J 2 - p 2 )$(p)]- 1/2 9[( J 2 - p 2 )Hp)]/dp, 



p= [-^{3} sing + cosg]^/±(J 2 - p 2 )<5>{p) (4.4a) 

and more simple in comparison with Eqs. (4.3) nonlinear Bloch-type equations obtained from 
Eqs. (4.3) by the substitution 

XjH cq -> XjH cmf = 

(2»{^}[^(ito)] 1/2 , -29f{»}[^(i^)] 1 / 2 , a + [»{5}yi - ^}lfe)][^(lto)]- 1/2 ^(ito)/^b) (4.46) 

Note that these latter Bloch-type equations are equivalent to those obtained in [10] in terms 
of variables Vi(t) =< Vi(t) > and solved in terms of hyperelliptic functions. 

So, Eqs. (4.1)-(4.2) and their specifications (4.3)-(4.4) yield a tool for examining quasi- 
classical dynamics of SL(2) clusters within subspaces £([/«]). However, they are not suitable 
for such an analysis at arbitrary initial conditions or for time-evolution of uncoupled (char- 
acterizing single subsystems) dynamical variables that is often necessary in practice. At the 
same time Eqs. (3.3) and (3.6) enable to obtain appropriate QAs 

U c H \t) = £ Sy^oY exp( ~'*y V) ) \[h];v)(v; [k}\ S Y (£ ) = 
£E«P( " ft ? 1;t,) )E£ S fv (J;g,r)S}M9,r)\[k];f)(f';M (4-5) 

[k]v>0 a />0/'>0 

of evolution operators Un(t) when eigenenergies E([li];v) are given by Eqs. (3.3a)-(3.4) or 
Eqs. (3.6). Evidently, L{ [/^-restrictions P^U^it) of such evolution operators (4.5) are 
distinguished from evolution operators U l H .^(t) associated with solutions of Eqs. (4.1a). 

Substitutions of Eqs. (3.3a)-(3.4) or Eqs. (3.6) in Eqs. (4.5) enable us to calculate 
appropriate QAs for time- dependences 

<F(t)>=Tr[U%(t)pU%(t)F] (4.6) 

of any dynamical variables F where p is a density operator for an initial quantum state. For 
example, inserting in Eq. (4.6) ordered exponentials of coupled (Yi) or uncoupled (original) 
dynamical variables one may get (after an appropriate Fourier transformation) formulas de- 
scribing dynamics of different (associated with GCS of both SL(2) and dynamic symmetry 
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groups of subsystems) types of Q-,P- and Wigner quasiprobabilty functions which are widely 
used for visualizing features of systems under study [14,26]. Note also that, due to Eqs. (2.5), 
the first line in Eq. (4.5) is more suitable for using Eq. (4.6) with F = F({Yj}) whereas 
the second one is more relevant for calculations with F depending on uncoupled dynamical 
variables. 

5 Applications to a class of quantum-optical models 

In this Section we manifest a physical meaning of general results above on ceveral concrete 
models which are widely applied in quantum optics, laser physics and quantum electronics 
[3,4,8,18,19,25]. Specifically, as was shown in [11], a natural area of applications of the sl p d(2) 
formalism is provided by quantum models with Hamiltonians 

Hjh = ]T Uia+ai + g'(at) m (a 2 ) n + g'* ( ai ) m {a+)\ n < m, (5.1a) 

m 

H 2 /h = ^uJidfdi + ujqOq a o + g'{ a i--- a m)( a o) n + g'*(ai...a m )(aQ ) n , n <m, (5.16) 
i=i 

N 

H 3 /h = u.afa, + 5>„(i)e/2 + g'a + (i)(a 1 ) n + g'*a_(i)(a+) n ] (5.1c) 
i=i 

where g' are coupling constants, Oj, af are boson operators describing field modes with fre- 
quencies Ui, <J a (i) are Pauli matrices, e is an energy difference of two level atoms and non- 
quadratic parts of Hi describe different multiphoton processes of scattering and frequency 
conversion (Eqs. (5.1a,b)) as well as the matter-radiation interactions in n-photon point-like 
Dicke models in rotating wave approximation (Eqs. (5.1c)). Note that in applications, one 
considers, as a rule, models (5.1) with n — 0, 1 that correspond, respectively, to semiclassical 
or completely quantum versions of models under study [3,4,8,18,23,25]. 

Appropriate Hilbert spaces L(Hi) are multimode Fock spaces Lp = Span{\{rii}) = 
Yli[ni\]^ 1 ^ 2 (a'l) ni \0 >} for models (5.1a)-(5.1b) whereas for models (5.1c) L(H 3 ) are direct 
products of single-mode Fock spaces Lp and "atom" spaces L a = Span{\j, fi; {j int } >)} 
where \j, [/,; {jmt} > are the basis vectors of irreducible representations of the "atom" group 
SU(2) a (with generators S a = J2iLi&a(i)) which are obtained from one- atom basis states 
|± > (i) with the help of the generalized Wigner coefficients and {ji n t} are sets of the SU(2) a 
intermediate angular momenta labeling basis vectors of the irreducible representations of the 
symmetric group Sn and being integrals of motion [11]. 

Hamiltonians (5.1) are expressed in the form (1.1)-(1.2) with the help of introducing 
sl p d{2) dynamic variables Vo,V + ,V- = (V+Y and integrals of motion Rj via a generalized 
Jordan- Schwinger mapping [10] given for Hi, H 2 , H 3 respectively as follows [11]: 

V = (a^ai — a^a 2 ), V+ = (a^) m (a 2 ) n , Ri = (na^ai +ma^a2), (5.2a) 

/// • // ' m + n 

j m 

V = — — (Y, a t a i - a t «o), V + = a+...a+ (a ) n , 
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1 1 m 

R k = ■ — (afa j -at 1 a j+1 ),k=l,...,m-l, R m = ■ — {n^af ai + ma^aQ), (5.26) 

m + n J J m + n ~[ 

1 N N n N 

^o-E^), V + = Y,° + (i)M a , i?i= gE^oW + afai (5.2c) 

i=l i=l i=l 

The structure polynomials ^(Vo) are determined with the help of Eqs. (5.2) (and defining 
relations for a(i), a + {i), a a (i)) from Eq. (1.3) which is valid for all L(ifj). Then for Hi, H 2 , H 3 
one finds, respectively, 

*(V ) = (mVb + Ri) (m) (Ri - raV + n) (n) , (5.3a) 

m—1 

*(V ) = [R m - nV + n]^([R m - (m + n) £ ii*]/™ + K,) ^ iV 2 . . . JV m _ ls 



i=i 

j m— 1 m—1 

N k = — [R m - (m + n) V iifc] + y + (m + n) V ifc, A; = 1, . . . , m - 1, (5.36) 

m i= i t=k 

*(V ) = [C 2 (2) - Vo (2) ][i?i - nV + n] (n) (5.3c) 

where C 2 (2) = E + E_ + (£ /2)( 2 ) is the Casimir operator of the "atom" su(2) algebra and 
[C 2 (2), V a \ = 0, = A(A - 1)...{A - B + 1). 

The subspaces £([/»]) in Eq. (1.4) are generated by the lowest vectors > which are 
given for different Hi as follows 

Hi: 11 

\[k] >= \{ni = re, n 2 = si), Z = (k — s), h — (nre + ms), 

m + n m + n 

Rj\[k] >=lj\[k] >, re = 0, l,...,m- 1, s = 0, 1,..., (5.4a) 

H 2 : 

m 

>= |{ni = Ki,n 2 = re 2 , • • • ,n m = K m ,n = s}), ]Jre; = 0, 

i 

lo = — - — (Vre, - s), l m = — 1 ~ — (nV/tj + ms), 4 = — — k = 1, . . .,m - 1, 
m + n J m + n ^ m + n 

^■IfZi] >=Z J -|[Z i ] >, re* = 0,1,.., s = 0,l,..., (5.46) 

I [h] >= \{ni = re}) |j, /i = -j; {j'^t} >, = ~h h = K- nj, 

Rj\[li] >= >, re = 0, 1, < j < N/2, (5.4c) 

where |{nj}) are standard Fock states and \j, —j; {jint} > are lowest vectors of irreducible 
representations of the "atom" group SU(2) a . From Eqs. (5.3)-(5.4) it follows that we have 
compact versions of algebras sl p d{2) in all cases except for models (5.1a,b) with n = 0. 

Eqs. (5.3)-(5.4) yield requisites for specifications of general results of Sections 3 and 4. 
However, for the sake of simplicity of our exposition, we restrict ourselves by considering 
certain simple examples which elucidate main features of new QAs and, simultaneously, will 
provide a base for further investigations of the most spread in applications models and phys- 
ically important cases (n < \,m < 3 in Eqs. (5.1a)-(5.1b) and n — 1 in Eq. (5.1c)). 
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Example 1. "Cluster" mean-field energy spectra in models (5.1) with n — l,m — 2,3. 
Inserting Eqs. (5.3) in compact (su(2)) versions of Eq. (2.1a) and using Eqs. (5.4) one finds 
for the polynomials $(lo) from Eqs. (2.1a) the following expressions 

Hl : = (mK °|;; J 1 + + 'yol' I)lm ' - J = sl2 > x>m) = MX (5 - 5a) 

m m m 

H 2 : <f> m (Y ) = U'(Y + J + Ki + 1) = n'(y + _ + Ki + l), n«i = 0, (5.56) 
i=i i=i " i=i 

H 3 : $(y ) = max( K) 2j) - J - F = { Zi'Jy^ > J = min 0'> f ) (5-5c) 

where it is also taken into account that d(M) = 2 J + 1 and simultaneously = s + 1 

for = 1,2 and d([li\) = min(2j, k) + 1 for H 3 ; besides, the numerator in Eq.(5.5a) 

always contains the factor (J + 1 + Y ) due to the definition of the symbolic powers 
and n^Li ' m Eq. (5.5b) means that in the product the term with = is omitted. Then, 
inserting Eqs. (5.5) in Eqs. (3. 3), (3. 4), (3. 6) and Eqs. (4.3)-(4.5) and using also Eqs. (5.4) one 
can obtain appropriate specifications of QAs above for energy spectra, Bloch-type dynamical 
equations and evolutions operators and to examine their features depending on characteristic 
parameters of models under study. 

However, postponing such a detailed analysis for further publications, we only write down 
appropriate specifications of Eqs. (3.6a) in the resonace cases (a = in Eq. (2.2)) and at 
m = 2 in Eqs. (5.5a,b) when Eqs. (3.6b) are solved analytically yielding 



. , 1 2k + 1 , n ,2/c + l w 2« + l ^ 
cos 2r = c( S , «) = - ^— + 1 - 2^1 + (— )( — + I! 

1 9k -+- 1 

w^(=^-^-l)(fors>2« + l) (5.6a) 



1 2k + 2 , / ,K+1 W K+1 

cos2r = c s, re) = - + 1 - 2W 1 + ( )( + 1) 

3 \ s V s s 

~I(^±l_l)( fors>K+ l) (5. 6& ) 

cos 2r = c(j, re) = J (l - 2//(re, j) + 2Jl-(j,(K,j) + (//(k, j)) 2 ) , //(re, j) = max ^y ) 
3 V v / mm(re, zjj 

(5.6c) 

for Hi, H 2 , H 3 respectively. Then, with the help of Eqs. (5.5)-(5.6) eigenenergies E cm f(\l^\; v) 

in the "cluster" mean-field approximation (3.5) are given as follows, 

Hi: 

E cmf {[h];v)/h-C = E cmf (K,s;v)/h- (k + 2s)ui = 



-\g'\(s - 2v) sin2r v /2((-s + 2v) cos2r + s + 2k + 1) 



-\g'\{s - 2v)J(l - [c(s, re)] 2 )2((-s + 2v)c(s, k) + s + 2k + 1), oo 2 = 2ui (5.7a) 
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H 2 : 

E cmf ([k};v)/h -C = E cmf (K, s; v)/h - km + k 2 uo 2 + s(wi + uo 2 ) = 
-\g'\(s - 2v) sin2ryj(~ + v)c(s, k) + & - + k + 1) = 

-\g'\(s - 2v) ] /(l - [c(s, kW)((- S - + v)c(s, k) + ^ + k + 1), 

K\K 2 — 0, /c = max(/ci, /c 2 ), c^o = ^i+^2 (5.76) 

E cm f(ik}; v)/h -C = E cm f(K,j; v)/h - (« - j>! = 
— |</| (s — 2f ) sin 2ry / max(/t, 2j) — J — (— J + t> ) cos2r = 

-|y'|(s - 2u)^/(l - [c(s, /t)] 2 )(max(/t,2j) - J - (-J + u)c(s,k)), e = wi, J = min(j, -) 

(5.7c) 

Evidently, Eqs. (5.7) manifest explicitly an essentially nonlinear dependence of energy 
levels E cm f([li\; v ) on their both s-u(2)-invariant (n,s,j) and non-invariant (v) labels unlike 
standard QAs obtained by means of using in Eqs. (3.1) GCS associated with dynamic sym- 
metry algebras of subsystems. Indeed, using in variational schemes (3.1)-(3.2) Glauber's CS 
Ilj T^{.Oii)\{ni}) , V{p.i) = exp(aiaf — a*cii) for models (5.1) and, additionally, "atomic" SU{2) a 
GCS exp(££ + — £*£_)|j, {ji nt } > for models (5.1c) as trial functions , one finds the following 
analogs of Eqs. (5.7) for such simplest QAs 

#i: E smf (n u n 2 )/h = ^(m + 2n 2 ) + A^i, |(/|), u 2 = %j x (5.8a) 

H 2 : E smf (iii, n 2 , n )/h = ui(ni + n ) + ^2(^2 + n ) + A. 2 (ui, oo 2 , \g'\), uj = co 2 + o;i (5.86) 

H 3 : E sm f(n 1: pi)/h = u x n x + fjn(u u \g'\) + A 3 (^i, \g'\), e = u 1 (5.8c) 

where Aj(. . .) are constant (for whole L(H)) energy shifts and f2 is an efficient frequency. 
Evidently, energy levels (5.8) depend linearly on Gj ni ,-noninvariant labels n« arranged on 
multidimensional lattices that provides multiperiodic dynamical regimes. Other ordinary 
QAs [4,8,25], e.g., obtained with the help of GCS of partially coupled subsystems, lead to 
similar results (as it is seen, in fact, from comparisons of Eqs. (5.8b) and (5.8c)). 

Example 2. "Cluster" mean-field energy spectra in models (5.1a,b) with n = 0, m = 3. In 
this case models under study yield so-called parametric approximations for models of the first 
example with m — 3. Besides, we have noncompact versions of sl p d(2) algebras because all 
subspaces L ([li]) are infinite-dimensional. This, in turn, causes an ambiguity of determining 
the parameter J in the generalized Holstein-Primakoff mappings (2.1) on subspaces £([/«]) 
because a polynomial character of $(Y ) = ^(Y + l — J + 1)/(J + Y )(—J + 1 + Y ) is 
provided by two values of J on each subspace £([/«]) that requires to add a choice procedure 
of J to Eqs. (3.1). However, we restrict ourselves by writing down analogs of Eqs. (5.5a,b), 

Hi : $ 3 (>o) = ^ ~^yI^+Y) = 27(F ° + Al(K ' Jk)) ' K = °' 2 ' (5 ' 9a) 
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H 2 - 



$ 3 (Y ) 



(-J+1 + Y )(J + Y ) 



Y + X 2 ({K i },J Ki ), 



H Ki = 0, (5.96) 
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which, nevertheless, manifest differences of parametric QAs from those given by Eqs. (5.5) 
due to linear and quadratic forms of $ 3 (Yo) in these cases (constants Aj(. . . , J...) are easily 
determined for chosen values J Ki ). A more detailed analysis of such comparisons will be 
given elsewhere. 

6 Conclusion 

So, we have obtained new approximations for energy spectra and evolution operators as 
well as nonlinear Bloch-type dynamic equations for models (1.1) (and (5.1)) by means of 
using the mapping (2.1) and standard variational schemes [17,5] with the SL{2) GCS as trial 
functions. They may be called as " cluster" (or correlated) QAs owing to taking into account 
strong quantum correlations between interacting subsystems. These approximations may be 
used to calculate in models of the (5.1) type time evolution of different quantum-statistical 
characteristics and quasidistributions (cf. [8,14]) and to find bifurcation sets and solutions 
of nonlinear Hamiltonian flows determined by Eqs. (4.1) and (4.5) (cf. [5,23]). In this way 
we hope to reveal in these models new cooperative phenomena and dynamical regimes (due 
to quantum correlations between subsystems) by analogy with those found in [9,18,22,23,27] 
and many other papers by using standard QAs; herewith different QAs above are expected 
to elucidate the role of such correlations depending on a choice of initial quantum states and 
paprameters of models under study (cf. [22,23]). 

However, from the practical point of view for this aim it is desirable to modify and 
to simplify Eqs. (3.3) and (4.5) by using different properties of the hypergeometric func- 
tions 2-P\(a, b] c; x), including their integral representations and asymptotic expansions [20,28]. 
(Specifically, in such a way one can express spectral functions Ej(r; J; v) in terms of the hy- 
pergeometric functions 4 F 3 (...; 1) which are proportional to the si (2) Racah coefficients [12].) 
Along this line it is also of importance to get estimations of accuracy of QAs obtained and of 
their efficiency in comparison with other approximations (e.g., considerd in [8,9,10,18]). One 
way to do such estimations is in comparisons of these QAs with appropriate computer cal- 
culations (cf. [18,27]) and another one is connected with using the "energy-trace" proximity 
meausures (3.7). 

Another line of further investigations concerns developments of mathematical aspects of 
the work. Indeed, results of Sections 3,4 correspond to picking out "smooth" si (2) factors 
<SV(£o) = ex P(£o^+ — m exact (generally, not "smooth") diagonalizing operators S 

determined by Eqs. (2.9), (2.12) and, when using Eqs. (4.1a), in evolution operators {/#(£) 
determined by exact evolution equations given in [11]; besides, Eqs. (4.5) yield another type 
of QA for evolution operators U H (t). All these QAs can be used as initial approximations 
in iterative schemes of constructing exact solutions which are similar to those developed to 
examine nonlinear problems of classical mechanics and optics [29] or in search of suitable 
multi-parametric improvements of variational schemes used by introducing "form-factors" 
with extra fitting parameters in original trial functions. It is also of interest to develop 
methods of obtaining simple formulas for exact solutions of tasks under consideration in 
order to compare with them results of approximations found above. At present one may to 
point out, at least, three promising ways along this line. 
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One of them is in simplifications of integral solutions obtained in [11] for both evolution 
and spectral tasks. The second way, leading to solving singular differential equations, is 
connected with using two conjugate differential realizations of sl pd {2) generators V a [10,12]: 

V + = z, V = zd/dz + l , V- = z' 1 ^ (zd/dz + l ), (6.1a) 

V- = d/dz, V = zd/dz + l , V+ = ^(zd/dz + ^(d/dz)' 1 (6.16) 

which are, in turn, related to realizations of sl p d(2) generators V a by quadratic forms in 
si (2) generators Y a taken in the coherent-state representations (cf. [30,15]). (In fact, these 
realizations were used implicitly when obtaining exact integral solutioms [11].) For exam- 
ple, when the structure polynomial ty(x) has the third degree (as, e.g., in models (5.1) with 
n = 1, m = 2), the realization (6.1b) reduces original tasks to solutions of the Riccati equations 
[12]. In this connection one may consider the hypergeometric functions 2 Fi(a,b;c;x) deter- 
mining QAs obtained as specific asymptotics of new classes of special functions determining 
exact solutions that opens a possibility to use the techniqe of asymptotic expansions [28] for 
finding latters. Last (but not least!) way is due to interrelationships between sl pd (2) algebras 
and certain g-deformed algebras mentioned in Section 2 that enables us to use for purposes 
formulated above techniques of g-deformed algebras and g-special functions, in particular, 
g-exponents defined with the help of the coherent states map of the paper [21]. Evidently, a 
progress in solving all these problems will promote to an extension of the orbit type GCS con- 
cept [16] and, simultaneously, to a more fine description of "classical" phase spaces associated 
with dynamic symmetry algebras sl p d(2) (cf. [31]). 
The work along these lines is now in progress. 
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